function [varargout] =Bmatic(varargin)
A=varargin{1};
AA=varargin{2};
p=size(A);
p=p(1,2);
for i = 1:p,
   s=A(1,i)/AA(1,i);
   c=A(2,i)/AA(2,i);
   AA(:,i)=AA(:,i)*s;
   ss(i)=s;
   cc(i)=c;
end
ss
cc
varargout{1}=AA;